Data_2002 = read.csv ("~/Desktop/PM 566/PM566-Labs/PM2.5_2002_Data.csv")
Data_2022 = read.csv ("~/Desktop/PM 566/PM566-Labs/PM2.5_2022_Data.csv")PM 566 Assignment 1
Assignment Details
In this assignment, I will be using Environmental Protection Agency (EPA) air pollution data to determine whether or not daily concentrations of PM2.5 have decreased in California from 2002 to 2022.
Step 1
Read CSV into Dataframe
The data for 2002 includes 15,976 observations (rows) of 22 variables (columns). The 2022 data has the same 22 variables (columns), but instead has 59,756 observations of each (rows).
Step 2
Combine 2002 and 2022 Data Into One Dataframe
Combined_Data <- rbind(Data_2002, Data_2022)Create New Year Column
Combined_Data$Date <- as.Date(Combined_Data$Date, format = "%m/%d/%Y")
Combined_Data$Year <- format(Combined_Data$Date, "%Y")Rename Key Variables
names(Combined_Data)[names(Combined_Data) == "Daily.Mean.PM2.5.Concentration"] <- "Daily_PM2.5"
names(Combined_Data)[names(Combined_Data) == "Daily.AQI.Value"] <- "Daily_AQI"Step 3
Map of Collection Sites
Although the monitoring sites are spread throughout California, they seem to be more concentrated along the coast, as well as in/around major cities (i.e., Los Angeles, San Francisco, San Jose, San Diego). Too, there are relatively very few sites in Southeast California (Eastern regions of San Bernardino, Riverside, and Imperial counties).
sum(is.na(Combined_Data$Year))[1] 0str(Combined_Data$Year) chr [1:75732] "2002" "2002" "2002" "2002" "2002" "2002" "2002" "2002" ...Combined_Data$Year <- as.numeric(as.character(Combined_Data$Year))
unique(Combined_Data$Year)[1] 2002 2022Combined_Data <- Combined_Data[!is.na(Combined_Data$Year),]
Sites <- (unique(Combined_Data[,c("Site.Latitude","Site.Longitude")]))
dim(Sites)[1] 202 2library(leaflet)
pal <- colorFactor(c("lightblue", 'darkred'), domain = unique(Combined_Data$Year))
leaflet(data = Combined_Data) |>
addProviderTiles('CartoDB.Positron') |>
addCircles(lat = ~Site.Latitude, lng = ~Site.Longitude,
opacity = 0.01, fillOpacity = 0.001, radius = 1, color = ~pal(Combined_Data$Year))Step 4
Based on some quick Google searches, most of these daily PM2.5 values seem plausible. Annual averages for California (specifically, Los Angeles) fall around 9 ug/m3, and daily averages can be as high as 35 ug/m3 for the same areas.
We see values much higher than this in our dataset (upwards of 50-69 ug/m3). Still, these values may still be okay as events like wildfires can drastically impact daily PM2.5 concentration averages (e.g., the 2018 Camp Fire in Sacramento lead to a daily PM2.5 concentration of 263 μg/m3, the highest ever recorded in California).
We also see a number of negative values with our daily PM2.5 concentrations. After some more Google searches, I learned that can occur because of two main circumstances: either there is some issue with a measuring instrument, or a measurement is taking place while the atmosphere is extremely clean (approaching 0μg/m3) and there is some level of measurement noise.
After a quick skim of the data, I’m leaning towards thinking that this data set’s negative values are due to the latter explanation, as the majority of them do not exceed -1.0μg/m3.
There do not seem to be any missing values for our variables of interest.
Step 5
Exploratory Graphs
library(ggplot2)State Level Trends
boxplot(Combined_Data$Daily_PM2.5 ~ Combined_Data$Year,
col = "lightblue",
pch = 20,
main = "Daily PM2.5 Concentrations by Year",
xlab = "Year",
ylab = "Daily PM2.5 Concentration (µg/m³)",
names = unique(Combined_Data$Year))
library(dplyr)
Attaching package: 'dplyr'The following objects are masked from 'package:stats':
filter, lagThe following objects are masked from 'package:base':
intersect, setdiff, setequal, unionmedian_results <- Combined_Data %>%
group_by(Year) %>%
summarize(median_PM2.5 = median(Daily_PM2.5, na.rm = TRUE))
print(median_results)# A tibble: 2 × 2
Year median_PM2.5
<dbl> <dbl>
1 2002 12
2 2022 6.8At first glance, the median value for the 2022 data (6.8 µg/m³) seems to be just lower than the same measure for 2002 (12.0 µg/m³). However, the IQR for the 2002 data is much greater than for the 2022 data, made clear by the differences in height of the pink rectangles representing this measure. Too, the max for the 2002 data was higher than for the 2022 data, and the opposite was true for the minimum. Finally, and although both years had a notable amount of outliers, the 2022 data had significantly more than the 2002 data. This difference could be due to an increase of severe weather events like wildfires that can provide extremely high PM2.5 values like those seen in the 2022 data.
County Level Trends
library(ggplot2)
Combined_Data[!is.na("County") & !is.na("Year")] |>
ggplot() +
geom_point(mapping = aes(x = Date, y = Daily_PM2.5), color = "darkred") +
facet_grid(County ~ Year, scales = "free")
When evaluating temporal differences in daily PM2.5 concentrations across counties, we can immediately notice distinct differences between the 2002 and 2022 data.
For instance, most of the data for 2022 is bolder, or more dense, than the 2002 data. This is because there are four times as many observations for 2022 than for 2002.
A handful of counties, Tehama, Modoc, Madera, and Glenn, do not have data for either one or both years, meaning that temporal differences in PM2.5 concentrations cannot be evaluated.
While the majority of PM2.5 concentrations maintain the same shape/trends between the two years, a number of counties’ concentrations look very distinct when compared across the two yearly data sets. More specifically, El Dorado, Mariposa, Mono, Nevada, and Placer counties have significantly higher values at various points throughout the 2022 year compared to the same time periods in 2002. These increases/peaks are likely due to increasingly-frequent extreme weather events like wildfires.
Site Level Trends
The site in Los Angeles that I chose to evaluate is the Los Angeles-North Main Street Station.
Combined_Data <- Combined_Data %>%
mutate(LA_Site = ifelse(Local.Site.Name == "Los Angeles-North Main Street", Local.Site.Name, "NA"))
LA_Site <- Combined_Data[!is.na(Combined_Data$LA_Site), ]
Combined_Data[!is.na("LA_Site")] |>
ggplot() +
geom_point(mapping = aes(x = Date, y = Daily_PM2.5, color= "LA_Site"), color = "lightblue") +
facet_wrap(~ "LA_Site", nrow = 1, scales = "free")+
labs(title = "Daily PM2.5 Concentration for Los Angeles-North Main Street Station, 2002 versus 2022")+
labs(x = expression("Date"), y = "Daily PM2.5 Concentration (µg/m³)")
The data for the LA North Main Street Site is reflective of the overall data for all California Sites. Although the 2002 data is more spread out than the 2022 data, the latter has significantly more outliers. For instance, the highest outlier for the 2002 data seems to be just over 100 µg/m³, whereas the max for the 2022 data is just above 300 µg/m³.